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q Abstract. Probability distributions defined on the half space are known to be quite different 
O 

from those in the full space. Here, a nonextensive entropic treatment is presented for the half 

•l-H 

^ space in an analytic and self-consistent way. In this development, the ordinary first moment of 
the random variable X is divergent in contrast to the case of the full space. A general v-th 
moment of X is considered as a constraint in the principle of maximum Tsallis entropy. The 
infinite divisibility of the distribution with an arbitrary v > and convergence of its /V-fold 
convolution to the exact Levy-stable distribution is discussed in detail. A feature of this 
derivation is that the Levy index is related to both the values of v and the index of 
nonextensivity. 
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It is known [1,2] that the properties of probability distributions define in the half space are 
quite different from those in the full space. Physically, such distributions are important in cases 
where several of the fundamental variables describing systems are bounded from below. Typical 
examples include time in any irreversible process and the Hamiltonian of a stable system. 

In recent years, the Levy distribution has been invoked to study many physical phenomena, 
such as anomalous diffusion [3], power-law responses of materials [4], aging in glassy systems 
[5], and so on. (See Ref. [6] for other topics.) The Levy-type distribution in the full space has 
been analyzed based on nonextensive statistics [7] with the Tsallis entropy [8]. The distribution 
derived therein has the same asymptotic behavior as the Levy distribution for a large value of the 
random variable. Combined with the Levy-Gnedenko generalized central limit theorem [1,2], 
this is expected to be a reliable procedure for obtaining the Levy distribution based on maximum 
entropy principle. Though the Levy distribution in the half space has been discussed in the 
literature and has been used in physical contexts [4,9], the consideration of this problem in the 
light of maximum entropy principle has not been given so far. The purpose of this paper is to 
provide such a demonstration for the half- space Levy-type distribution to fill this gap. 

The Levy-stable distribution in the half space (0, °°) is defined by [2] 

L r (x) = ^-1 dk XL (k)e- k \ (1) 
2k j 



In this equation, % L (k) is the characteristic function, which is given by 



X L (k) = exp j-<a| k \ 1 exp 



ie(k) 



On 



(2) 



where a is a positive constant, y the Levy index, 6 a constant satisfying 1 6 \ < y and 
e(k) = k/ 1 k | the sign function of k . In contrast to the Levy distribution in the full space whose 
index is in the range < y < 2 , in the half space the index is confined to 
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0< 7<1. 



(3) 



This distribution has the following asymptotic behavior for large values of x : 

L r {x)~x- l -\ (4) 

and therefore all of the ordinary positive-integer moments, (X n ^=^ dx x" L y (x) 

(n = 1, 2, 3, L ), are divergent. This should be compared with the case in the full space where the 

(only) ordinary first moment can be finite. (In the notation used here, x is the value of the 
random variable X .) It is known that the Levy distributions |L r (^:):0< 7<l} form a stable 

class in the sense that the sum of N scaled independent random variables {* ; 2 L N 



x = — — '*" 2 ~*~^ ~^~ Xn (5) 
B N 



obeys the same distribution as Levy's. The scaling factor B N has to be chosen in such a way that 
the limit distribution is independent of the number Af of convolutions. Suppose each x i obey the 
Levy distribution in eq. (1). Then, the distribution of x is given by 

Lf(x) = B N (L r *L *L r ){B N x), (6) 



where 

x 

(/**)(*) = Jdfc' /(*-*') *(*')• (7) 

o 



The characteristic function of liy\x) is calculated to be 

X [ N) (k) = ] dxe^L^(x) 
o 
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f 



k 



= X L 



(8) 
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In calculating the inverse Fourier transformation, note that L y (x) = for x < . Thus, we find 



B N = N 



(9) 



where, without loss of generality, we may set the constant of proportionality to be unity. 

The derivation of the Levy distribution in the full space based on the ordinary maximum 

entropy principle with the Boltzmann-Shannon entropy turned out to be unsatisfactory as shown 

in Ref. [10]. This difficulty was overcome by using Tsallis' nonextensive generalization of the 

maximum entropy principle [7] with the constraint on the generalized second moment defined by 
(X 2 ^ =| dx x 2 [p(x)] g 1 1 dx [p(x)] q , where p(x) is the distribution to be determined by this 

principle and q is the index of nonextensivity. So far, this formulation is given only in the full 
space. In what follows, we examine the corresponding problem in the half space. In this 
consideration, to keep the discussion general, we treat a class of the v-th moments in place of 
the generalized second moment as a constraint. 

The nonextensive statistics formalism is based on the Tsallis entropy [7,8] 




(10) 



where a is a scale parameter and the distribution p(x) is normalized to unity 




(ID 
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We maximize S q [p] under the constraints on the normalization in eq. (11) and on the 
generalized V-th moment defind as 
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I dx x v \p{xj\ ' 



X v ) =± 



\ dx[p(x)Y 



(12) 



where v is an arbitrary positive number. The resulting distribution function is found to be 



p(x) 



(13) 



where Z q (f3) is the normalization constant, = f3 1 o q 1 with f3 the Lagrange multiplier 
associated with the constraint in eq. (12) and c = | dx [p(x)Y . In the present work, since we 
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are interested in the Levy problem, we focus on the distribution with a long tail. Then, q should 
be larger than unity. The normalizability condition requires q < v + 1. In addition, the divergence 
of the ordinary first moment (X), places another condition q > v/2 + 1. Therefore, the range of 

interest is 



+ l<q< v + 1. 



(14) 



It should be noted that this condition guarantees the finiteness of (X v ) . The quantities c and 

\ / q " 

(X y ^ have to be calculated by using p(x) self-consistently. The normalization condition on 
p{x) leads to the identical relation 



= [z,(«] 



\-q 



(15) 



Let us rewrite eq. (13) in the form 



p(x) = 



Z q (P) 



( \ 
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(16) 



where 



s = 



q-l 



(17) 



(18) 



The range of q in eq. (14) is now written in terms of s as 



1< vs<2. 



(19) 



From eqs. (16) and (18), consistency requires that £ v be positive. That this turns out to be 
consistent with the above inequality will be ascertained subsequently. Using the distribution 
function in eq. (16), we have the following results: 



B 



n i 

—,s 

V V V 



(20) 



X v ) = 



( sc A 



q vZ,(j8) 



B 



1 



+ 1, s 

V v v 



(21) 



where B(w, z) is the beta function. Combining these results, we obtain 



X v ) 



(22) 



Substituting this into eq. (18) and taking eq. (19) into account, we confirm that £ v is positive as 
promised above. Finally, the distribution is found to be 



6 



B 



V V V 



(23) 



Now, we are in a position to examine the properties of convergence of this distribution to the 
Levy- stable distribution after iV-fold convolution for large N. The characteristic function of the 
distribution of x in eq. (5) is 



X m (k) = ]dxe ikx B N (p*p*L *p){B N x) 



X 



\ B Nj 



(24) 



where x(k) is the characteristic function of p(x). Following Ref. [1], we express %{kl B N ) as 
follows: 



X 



l + ^dx 



( ■ 



exp 



ie(k)\ k \ x 

B N j 



-1 



p(x) 



i / i i \ Vi-l 



B 



fl 1 



B A 



,5 — |v-jvy o 



]dy 



exp[ie(k)y ]-l 



+ - 



B 



fi i , .. 

,s- |v-wy 

V V V 



x J dy {exp[z'f(fc)y ] - 1} 





l 




1 






V 


» y vs 


y v + 






V B N J 







(25) 
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To evaluate the integral, we use the following result given in Ref. [1] (p. 169): 



e ±iy - 1 



f 



,1+a 



an 



-exp mi ] | L{a) \, 

V 2 



(26) 



where L(a) is given by 



L (a) = ]^(e- y -l)<0 (0<a<l). 

0? 



(27) 



Equation (26) is then calculated to be 



f 7 A 



X 



\ B NJ 



I- a 



V B NJ 



exp 



le(k) 



On 



+ c 



N ' 



(28) 



where 



v^-'|L(v,-l)| ^ 
a = 7-r — 1 > 0, 



B 



1 1 

V. v v 



(29) 



= 1- vs. 



(30) 



The term C N in eq. (28) corresponds to the thrid term in eq. (25) and is the higher-order 

correction, which vanishes faster than the second term in the limit N — > °° . Taking the logarithm 
of X (N) (k) m e£ l- (24) and comparing it with eq. (2), we have in the limit N — > °° 



ln;r w (fc) = iVln* 



= -a I k I r exp 



ie(k) 



0n~ 



(3D 
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where the Levy index y is identified to be 

y=vs-l, (32) 

and B N is the same as that given in eq. (9). From eq. (19), equation (32) is seen to be consistent 
with eq. (3). Equation (31) shows that the nonextensive statistical distribution in the half space is 
infinitely divisible for an arbitrary v > and approaches the Levy-stable distribution by the N- 
fold convolution for large N. Equation (32) reveals an important feature of the approach 
presented here in that the Levy index is related to both the order v of the moment and the index 
of nonextensivity q = 1 + 1/ s. 

In conclusion, we have examined the nonextensive statistical formalism to investigate the 
Levy distribution in the half space. We have considered the generalized V-th moment of the 
basic random variable and have found the self-consistent solution of the nonextensive statistical 
distribution. We have shown that this distribution is infinitely divisible for an arbitrary v > 
and, in full conformity with the generalized central limit theorem, the iV-fold convolution of the 
distribution converges to the Levy-stable distribution in the large-A^ limit. 
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